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ABSTRACT 


The barotropic vorticity equation is utilized to investizate 
the steady topographic disturbance of an initially uniform zonal 
current. Solutions are obtained on the beta-plane for various dis- 
continuous bathymetric profiles, It is shown that a smooth slope 
does not appreciably modify the results, Stability criteria are 
derived that define the values of the flow parameters for which 
the streamlines become discontinuous, 

The solutions are evaluated for several observed disturbances 
of oceanic currents. The velocities implied by the observed wave- 
lengths are physically realistic, but the discrepancy between 
observed and calculated amplitudes suggests that the vertical 
structure of the velocity profile and the baroclinicity of the 


ocean are important. 





CHAPTER I 


INTRODUCTION 















In most studies of the large scale circulation patterns of the 
oceans, the sea floor is assumed to be either flat or infinitely deep. 

Such simplifications are made to facilitate the investigation of certain 
classes of flows in which it is suspected that topographic effects are 
nonexistent or of minor importance. That currents in the real ocean are 
to some extent topographically affected nas been pointed out by many 
authors. Perhaps the most universally cited example of such an effect is 
the wave-like structure in the dynamic topography of the Antarctic 
Circumpolar Current (ACC) where it crossesseach of several ridges (see 
‘Figure 1). 

Many contributions to the understanding of oceanic currents have been 
achieved through the application of the techniques and principles of 
dynamic meteorology. That —_ apparent similarities exist between oceanic 
and atmospheric flows was graphically demonstrated in 1951 when Fuglister 
and Worthington published the results of an intensive multiple ship survey 
of a Gulf Stream meander. Their illustration of the motion and incipient 
occlusion of the meander bears a striking resemblance to an earlier 
illustration (Berggren, et al., 1949) of the unstable waves observed at 

the 500 mb level in the atmosphere. It will therefore be found instructive 
to consider briefly the results of previous investigations of topographic 


effects on both oceanic and atmospheric flows. 





I-1. Oceanographic Investigations 
Elman (1923), in an attempt to account for the deflection of the 


Gulf Stream over the Southeast Newfoundland Ridge, considered the effect 

of bathymetric relief on a homogeneous current under a balance of pressure 

gradient, bottom aeaaitilon, and constant Coriolis forces. His results show 

the current deflection to be dependent on the slope of the bottom, but not 

the total depth, and give a net lateral displacement of the streamlines (to 
the right for flow over a ridge in the Northern Hemispnere). For flow over 
Successive ridges, Ekman's streamlines are 90 degrees out of phase with the 
topography, the maximum deflection occurring at the point of inflection of 

the bottom profile. 

In 1932 Ekman reconsidered the problem and included accelerations 
in addition to bottom friction, constant Coriolis force, and pressure 
gradients. Subject to the requirements that both the bottom slope and the 
variation in depth be small, he obtained a solution that gives streamlines 
in phase with the bathymetry, the maximum displacement to the right occurring 
over the point of minimum depth. 

In an attempt to rigorously investigate the effect of bottom friction 
on a topographically disturbed homogeneous current, Girtler (1941) utilized 
Prandtl's boundary layer theory. He simplified tne problem by considering 
a constant Coriolis parameter and a variation in depth that is small 
compared to the total depth of the fluid. His solution to the linearized 
equations reduces to that of Ekman (1932) provided that Gtrtler's friction 
parameters, including the boundary layer thickness, are properly 
identified with Ekman's "depth of frictional resistance." G&rtler extended 
Ekman's solution by including the possibility of the upstream flow 


approaching the topographic feature at other than normal incidence. He also 





; evaluated the frictional effects for several pnysically reasonable values 
of tne flow parameters and thereby demonstrated that such effects are 

' negligible for most atmospheric and oceanic currents. Gé&rtler pointed out 
though, as did Ekman (1932), that topographic effects in which inertia 
plays an important role may not properly be treated as independent of 
planetary effects when the scale of motion is such that the latter must be 
considered. This point is central to the investigation at hand. 

Sverdrup (1941) took a new approach to the problem of topographic 
effects. He neglected friction and considered that the effect of topographic 
features was to disturb the internal distribution of mass. Being unable to 
treat the problem rigorously, he nevertheless was able to demonstrate 
qualitatively that a bathymetric feature which rises from a layer of no 
motion into an internal solenoidal field and its associated Giprene will 
disturb that current in such a way that the physical sea surface will be 
depressed above the feature. This effect causes the current to deviate to 
the right (in the Northern Hemisphere), the maximum deviation being over 
the point of minimum depth, and then to return to its original configuration 
downstream. That the ain is true for flow over a depression is not 
clear, as there is no mechanism for disturbing the solenoidal field. 
Sverdrup's conclusions agree qualitatively with the observed features of 
the circulation in the Southern Ocean. 

Shtokman (1947), in a discussion of his subsequently published (1948b) 
theory of topographic effects on ocean currents, demonstrated an apparent 
serious discrepancy between the theories of Ekman, G8rtler, and Sverdrup 
and observed flow features of the South Atlantic and Caspian Sea. Considering 
Defant's (1941) chart of the "absolute" topography of the surface in the 


Atlantic, and relating this chart to the bathymetry of the South Atlantic 





between 25° and 50° South, Shtokman observed that the eastward currents 
are deflected poleward over shoaling water and equatorward over deepening 
water. Tnis is contrary to the observed deflection of the ACC. Since 
the direction of deflection in the earlier theories depended directly on 
the sign of the Coriolis parameter, there was no mechanism available to 
explain opposite deflections in the same hemisphere. 

Noting the above apparent discrepancy between theory and observation, 
Shtokman inferred that the Coriolis parameter, if considered constant, 
Should not be the primary consideration. On this basis, he adapted an 
earlier (1948a) derivation of the equations relating wind stress to lateral 
friction to the problem of a variable depth. For zero wind stress, the 
Signs of the terms in his equations give qualitative agreement with 
Sverdrup's (1941) conclusions. Shtokman's approximation of the complete 
equations including wind stress implies that the direction of the current 
deflection over a topographic feature depends on the sign of the local 
curl of the wind stress and not on the hemisphere. Thus, for shoaling 
water, a negative wind stress curl yields a deflection to the right, and 
vice versa. 

As Shtokman pointed out, his results agree with the observations he 
cited. Thus the absence of perturbations on Defant's 191 topography of 
the surface of the South Atlantic between 10° and 50° South is explained 
by the existence of maximum wind stress, and hence zero wind stress curl, 
in these latitudes. A later presentation by Defant (1961) of the dynamic 
topography of the same area indicates substantially the same structure 
except for the addition of an anticyclonic eddy at approximately l5°S - 
20°F (see Figure 2). If this is indeed an area of zero wind stress curl, 


then Shtokman's theory does not account for this eddy. 





The writer believes that a more fundamental objection to Shtokman's 
theory may be made on the basis of his method of introduction of a variable 
depth. The depth which Shtokman ultimately considers as representing the 
actual bathymetry was originally defined as that depth at which the 
transport components vanish. Charney (1955) and Fofonoff (1962) have 
shown that this depth of no motion is a function of the transport itself, 
and hence is of a nature that is fundamentally different from the 
bathymetry. 

Shtokman's refutation of the theories of Ekman (1923), Gertler (1941), 
and Sverdrup (1941) based upon consideration of the near shore currents in 
the Caspian Sea seems inappropriate in that those earlier theories specifi- 
cally excluded changes of depth that amount to a significant fraction of 
: the total depth. The earlier theories also excluded flows initially 
parallel to the isobaths, such as are considered by Shtokman (197, Fig. 5). 

Neumann (1960) attempted to combine analytically the results of 
Ekman (1923) and Sverdrup (1941). Neumann initially considered the steady 
flow in a baroclinic ocean to be maintained by a balance of Coriolis, 
friction (both wind stress and "internal" friction), and pressure gradient 
forces. He pointed out, as have Defant (1941) and others, that for 
homogeneous water, the variation of Coriolis parameter with latitude requires 
that all purely geostrophic currents be zonal, unless the depth increases 
from the equator toward the poles in proportion to the sine of the latitude. 
Neumann obtained approximate solutions for several classes of flow, including 
that of homogeneous water flowing over a ridge under the influence of a 
variable Coriolis parameter. This latter solution gives a Sverdrup- type 
deflection of the streamlines without the consideration of a solenoidal 


field. Neumann's solutions are based upon the requirement that the change 





of depth be small compared to the total depth, but the numerical examples 
given do not clearly meet that criterion. : 
Saint-Guily (1962) has extended his earlier (1959) treatment of the 
generalized Ekman problem to the case of a variable deptn. Subject to the 
assumptions of vanishing friction and moderate changes in depth, Saint- 
Guily obtained solutions for the steady zonal flow of homogeneous water 
over a sill that correspond to those of Neumann (1960) for a variable 
Coriolis parameter. Saint-Guily's initial disregard of the field 
acceleration terms in the equations of motion restricts the validity of 


his solutions to regions where the ratio of the relative vorticity to the 


Coriolis parameter (Rossby Number) is small. 


I-2. Meteorological Investigations 

AS was noted above, topographic effects on ocean currents have 
parallels in meteorological problems. Queney (1949) summarized previous 
theoretical treatments of topographic effects, including the results of 
nis own extensive treatment of the problem (19l)7) utilizing perturbation 
techniques. Queney considered the stability associated with the density 
stratification of an infinitely deep atmosphere, and it is therefore 
difficult to apply his results directly to a homogeneous ocean. Charney 
and Eliassen (1949) concluded that the real gravitationally stable atmosphere 
may be treated as an "equivalent barotropic atmosphere" whose properties 
coincide with those of the real atmosphere at a height of zero horizontal 
divergence (approximately 500-700 mb). Charney and kliassen were able to 
substantially improve their numerical prediction techniques by including 
a rough model of the effects of continental topography. Tne importance of 
such continental effects in the study of large scale atmospheric flows has 


been discussed by Bolin (1950). He noted the regular perturbations of the 


t nN 














pressure surfaces that occur on monthly (and longer) mean weather charts of 
both the Northern and Southern Hemispheres and argued in favor of a topo- 
graphic mechanism for the generation of such perturbations as opposed to a 
thermodynamic mechanism. Bolin linearized the equations and obtained 
approximate solutions for flow over an isolated mountain. 

Oi (1956) considered the problem of the frictionless barotropic flow 
of air over a mountain range under the influence of a constant Coriolis 
parameter. He was able to include accelerations in a manner similar to 
‘Bh (1932) and, again for small changes in the total depth of the fluid 
column, obtained approximate solutions showing the effect of a mountain 


range on an initially zonal flow. 


I-3. Initial Evaluation of the Problem 
All of the above investigations, with the exception of that of 
Shtokman, have essentially linearized the problem by considering only small 
changes in relief relative to the total depth of the fluid column. This 
may be shown (Bolin, 1950) to be equivalent to assuming that the relative 
vorticity of the fluid is smal] compared to the local value of the Coriolis 
parameter. In the ocean, with a mean depth of about 4 km, and in the 
equivalent homogeneous atmosphere, with a height of 8 km, it is apparent 
that there are many areas where such a linearization will be inappropriate. 
At low latitudes the relative vorticity may certainly approach or exceed 
. the value of the Coriolis parameter. In the oceans there are innumerable 
submarine features, covering a spectrum of size scales, that rise well above 
| the mean surrounding depths. These observations on the inherent limitations 
- of the perturbations techniques, coupled with the well known importance of 
the variation of the Coriolis parameter with latitude, provide the basis 


for the present investigation. 























In subsequent chapters, Gtrtler's (1941) conclusions will be 
accepted and it will be assumed that the flow of the barotropic ocean is 
frictionless, although it is recognized that friction must be important 

in some limiting sense. Indeed, G6rtler pointed out that frictional effects 
nay become important when the horizontal scale of the topograpnic feature 
under consideration is very much larger than the inertial radius of the 
flow. As the objective of this investigation is to study the combined 
effects of the variation of Coriolis parameter and finite changes in depth 
over a range of scales, subsequent assumptions will not restrict these 
essential features of the problem. 

The cencept of vorticity has become fundamental to the investigation 
of both atmospheric and oceanic flows. In its general form, the vorticity 
equation is a concise statement of one of the laws of fluid motion. For 
the class of flow under consideration its form facilitates investigation 
7 either an Eulerian or Lagrangian viewpoint. This property of the 
vorticity equation will be utilized in succeeding chapters. 

The deep currents in the oceans will be treated as though they were 
barotropic, not because this mode is necessarily strictly representative of 
in situ deep currents, but because the barotropic equations of motion are 
Bepler to treat exactly than are the baroclinic equations. The influence 
on this barotropic flow cf topographic features with the scales of features 
Bee tine in the real ocean will be investigated. 

Topographic effects observable in tne baroclinic flow (see Figure 1) 


must be related to those occurring in the barotropic flow, but this relation 


has not yet been determined. 





CHAPTER II 


FORMULATICN OF THE PROBLEM 


II-]. Basic Equations 


We consider the horizontal components of the equation of motion of 










mw -».-22 De, gym - LOE 
mo) Boe? "ae a 
together with the continuity equation 
JU, ove Ow 
yx * Dy Dz? (2) 


where u, v, and w are the velocity components along the x (eastward), 


y (poleward), and 2 (upward) axes. The Coriolis parameter (equal to twice 


P is the local pressure. The local acceleration of gravity is g and f is 


a 


the fluid density, both of which are assumed constant. 


mii-2. The Vorticity Equation 


The two components of equation (1) are cross-differentiated with 
respect to y and X respectively, and then subtracted to give, after 


rearrangement of terms, 


DF+ f) | u_ dvi Jw du _ dw dv 
ht Grn ES ee ee (oF 
where F = ov ~ = ,» the vertical component of relative vorticity. 


Equation (3) may be considered as the primitive vorticity equation, 


where the change in the absolute vertical vorticity of a parcel of fluid is 





10 











seen to be due to the local two-dimensional velocity divergence and the 
conversion of horizontal vorticity to vertical vorticity. 

If we now consider u and v to be independent of z, the vortex tube 
term in equation (3) drops out, and we may integrate equation (2) from the 
bottom of a column to the free surface to obtain 


Ju sue. 1 D(h +) 
Dooery “oad |) ee , 


where >| is the elevation of the free surface above its undisturbed position, 
and h is the height of the undisturbed free surface above the rigid bottom, 


or ground. Substitution of the above equation into equation (3) then yields 


me (if «ee 

& EE] -o. eI 
the equation cemmonly interpreted as expressing the conservation of potential 
vorticity along a trajectory. It is of interest to note that no assumption 


of steady state has been made in the preceding developnient. 


II-3. The Transport Streamfunction 


To facilitate the introduction of the volume transport streamfunction, 


the vertically integrated continuity equation may be written as 


[ + ls + n + ne 
dur h) + ge mw ar Ms - aa. 


For stationary conditions this equation is satisfied identically by a volume 


transport streamfunction defined such that 


at = _ u(r a) , and #¥ a vy +h); where Y = W(x,y). (33 


Henceforth, we will neglect the elevation of the free surface a compared to 


the depth h, both as a readily justifiable approximation and for expediency. 








PL 















In terms of the streamfunction y , the first integral of the 


icity equation, equation (4), is then 


he 
VC = + la constant (along a streamline). (6) 


955) and Fofonoff (1962) have shown that there must be a 

ional dependence of the Coriolis parameter and the depth h on the 
samfunction ¥'. This may be readily seen by considering that the above 
ation is applicable along an arbitrary streamline, and that the choice 
| specific streamline then determines the geographic path followed by a 
d column and hence its future depth and latitude. Substantial 
alytical difficulties arise when an attempt is made to express both the 


-hymetric relief and the variable Coriolis parameter as functions of the 


ream flow pattern. 
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CHAPTER II] 


SOLUTION FOR ZONAL FLOWS 


The difficulties mentioned at the end of Chapter II may be alleviated 













by restricting our consideration to flows that are strictly zonal upstream 
of the region of interest, and to topography that is independent of y. 


III-1. Simplification of the Problem 


o 


Specifically, we now establish the following conditions: at some 
point upstream of the region of interest, say at x = xX, the flow is 

defined by = - UHy, where H is a constant depth and U is the transport 
velocity (positive eastward); in the first region of interest, say x7Xp, 

the depth is a function of x only, h(x). In all regions the Coriolis para- 
meter is proportional to y, the coefficient of proportionality being 4, and 
with the origin of the y axis at the equator, f = “y. Denoting the upstream 
region by the subscript 1, we then have Y = -UHy,, hence f)(y,) = 


By (¥ 7) - - PY, /UH. Substituting into equation (6), 


1 
UH h(x) 


where YW * al (because equation (6) applies along a streamline). 


Rearranging and dropping the now superfluous subscript 1, 





h(x) Vy - 5 PEG) y - ~ n@(x)By . (8) 
















III-2. Matching Conditions 

The flow defined by equation (3) must satisfy continuity requirements, 
that is, it must match the upstream transport at x = x,. As the upstream 
transport is geostrophic, we might meet these requirements by matching 


meridional pressure gradients, but this is found equivalent to directly 


matching total transport, or transport components. The matching conditions 


are then 
oY . df - 6 
ox OX 
at x= xX, . ie 
d¥- a. - i 
s y o) o 


III-3. The General Solution 
Tne vorticity equation (3) may readily be reduced to an ordinary 
differential equation by utilizing the separation W(x,y) * y Z(x). 


Upon substitution and cancellation of y, equation (8) becomes 





2 
nox (x) = n(x) + EE ge - nhs (10) 
where pe = Ary ; 
The matching conditions (9) become 
J(x,) e ~ UH (11) 
Zix)) = 0 


To facilitate the solution of (10) subject to (11) we make the 


following transformation of the independent variable x: 
x 
Let (x) = G(s), where s » Bf nea, and 8, *= s(x, ) = 0, 


and define D(s) by: h(x) ® - UH“D(s). 





ly 










Jpon substitution, equation (10) becomes 


2 * Go = 1/Da(s) (12) 


- 
~~ fe 
oe ae 


tne matching conditions (12) becomes 


Go(0) = G\(0) = - UH 

| (13) 
G,(0) . G, (0) pe 
where the subscript 2 indicates the second region considered, that is, the 
region immediately downstream of the uniform zonal flow. Anticipating that 
we may desire to consider more than one downstream region, we may 
immediately generalize equations (12) and (13) to: 


G! + Gy * 1/D,(s) (1) 


G,(3,,) ’ Gn fen-1(™) 


me 2, 3, Lb, <a 
oh, fs, (a), 


t 
GC e I= 


no 


a 
where s(a) # ef h(x)dx 
O 


The complete solution of equation (14) is 


n 
ie) = 4, cos s, + B, sin 5, f sin(s, - t) dt (16) 
Sno | Dnt t) 


where the coefficients A, and B, are to be determined from equations (13) 
and (15). 

The solution (16) may be evaluated (if only numerically) for any 
bathymetric profile that is independent of y. It represents the topographic 


disturbance of a zonal current subject to the following assumptions: 


7 











Horizontal velocities are invariant with depth. 

The influence of friction is negligible. 

The convergence of the meridians may be neglected, and the 
meridional variation of the Coriolis parameter may be 
approximated as a linear function of distance from the equator 
(8-plane approximation). 

Displacements of the free surface from a level surface are small 
compared to the total height of the fluid column. 

The initial flow is zonal with zero relative vorticity. 


The flow, and hence the pressure field, is stationary. 
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CHAPTSR IV 


gVALUATION OF TH SOLUTION FOR SOME SIMPLE PROFILES 











when solution (16) is transformed back into Z(x), it is apparent 
that eastward flow over a varying depth is periodic, while westward flow 
S exponential, as p is real or complex depending on tne sign of U. 


Throughout this chapter, only eastward flow will be considered. 


Iv-1. The Bottom Profile is a Simple Step 

We consider first a simple step-shaped feature in the topography. 
For x<0O, h, = H; for x70, ho = constant. Then, transforming (13) and 
(16) back in terms of x, we have s * ropx, where ry * ho/H. and Ap 2 
-UH, By = 0. These values of the parameters give 


W(x,y) = y A(x) = - UY (1 - (1 - rp)cos ropx) . (17) 
rp 


It is apparent that #(x) will vanish at some point for rn%2. We 
may describe this situation as unstable because at such a value of x there 
Will be no zonal flow. The characteristics of solution (17) are best 
illustrated by reducing all streamlines to the single path given by 

. UY . - 

» (I= (I -"r5)cos ropx) 

The left hand member of equation (13) is shown plotted against ropx 
(Figure 3), for the stable case ro = 0.5 (Figure 3a), and for the unstable 
case ro = 2.5 (Figure 3b). In the unstable case the dashed streamline has 
been included to indicate the cellular structure extending to infinity at 


periodic values of x. This point will be further discussed in a subsequent 
chapter. 





af 

















It is instructive to make an estimate of the wavelength of the 
specific stable flow shown in Figure 3a. Taking o =» 2x 10723 em7+sec7* ; 
and arbitrarily choosing 1 cm/sec as a characteristic velocity of the deep 
oceanic circulation (only the order of magnitude of the velocity is impor- 
tant, the wave-length depending upon its square root), the wave-length is 
on the order of 280 km. For an atmospheric velocity of 20 m/sec, the wave- 
length would be on the order of 12,000 km, that is, of continental 


ai mensions. 


Iv-2. The Bottom Profile Resembles a Plateau 
We next consider a more complex example, that of flow over a plateau 
(or over a flat bottomed, vertical sided valley). We establish the 
conditions that for x<0, h = H as before, and for O0<x<a, h*= hp. For 
x7a, h3 = H, although the treatment is as easy with the final depth 
different from the initial depth. 
From the previous example, 

4, = By = B, * 03; Ap = - UH, 
and from the above definition of the topographic feature, 
mF F 


3715 ro * ho/H . os 


The boundary conditions on P(x) are then 


$3(a) * Bola) = - gs (1 - (1 - rp)cos ropa) , 


v 


; Cl 
23(a) a P>(a) = - Up (1 - ro) sin ropa . 


From equation (16), transforming back in terms of x, 


$3(x) = (A3 + UH)cos p(x ~ a) - UH + Bysin p(x - a) . (20) 
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Solving equations (19) and (20) for A3 and B3 and substituting back into 

equation (20) we obtain 

¥3(x,y) =y $3(x) = — UHy l(a/r, - 1)(1 - cos rgpa)cos p(x - a) a 
+ (1 - rg) sin ropa sin p(x - a) + 1 

Equation (21) describes the flow for all x>a and for ro either greater 

than or less than unity. 

In our analysis of equation (21) we will only consider values of the 

parameters ro and a that give stable flow in region 2. It does not seem 

realistic to match a solution that is unstable with a subsequent downstream 

flow. 

Proceeding now to the development of stability criteria for the 

downstream flow, by our earlier definition instability occurs for some 

value of x when the bracketed term in equation (21) vanishes, or: 

(1/ro - 1)(1 - cos ropa) cos p(x - a) + (1 - ro) sin ropa sin p(x - a) = -1. 

If we now let the coefficients of cos p(x - a) and sin p(x - a) be 

represented by C,(ro,a) and Co(ro,a) respectively, and divide the above 

equation by (cf + c2)2 we obtain 


Cy C ~1 


cos p(x - a) + e sin p(x - a) = ‘ 
(c2 oF : (+ eye (Fe eye 


By defining the angle g such that tan ¢ = C,/C,, and utilizing a simple 


trigonometric identity, this equation can be written 


o 
(ce + ce sin (g¢ + px - pa) = -l. 


Noting that ¢ and pa are phase angles, this equation will be satisfied at 


Some value of x only if ce 


2 
1* Cy as 








ie 


In terms of ry anda, the criterion for the instability of the 
















downstream flow is 


[ 2 
2 i: af 
(1 - cos rgpa)© + rj sin*rjpa 2 ers : (22) 


No further general simplification of the criterion is possible, equation 
(22) being transcendental in roe However it would appear that there is no 
value of a for which an ro cannot be found that will satisfy the above 
relation. As must be true, there is no instability for ro * l. A rough 
calculation indicates stability for ro between O.67 and 2. As both ro and 

a are topographically determined, the properties of the solution are best 
illustrated by example. 

We first investigate the flow when ro = 0.5, previously found to be 
Stable throughout region 2. The solution (21) for the downstream region may 
then be written 


ey . f + (1 - cos 0.5 pa) cos p(x - a) + 0.5 sin 0.5 pa sin p(x - ” pa 
B ) 


af 


Four representative breadths ‘a will be chosen. The first, 0.5 pa ® 7/3 
(or 77/3, etc.) locates the down-step in a region of S.&. flow; the second, 
at 0.5 pa = 5/3 locates the step in a region of N.E. flow. The third 
Peeadtit is chosen such that the step is at the point of maximum excursion, 
0.5 pa =7T; and the fourth at the point where the flow has regained its 
original latitude, 0.5 pa = 277, 

Figure l shows the normalized streamline for the first case. The 
flow over the plateau is the same as that shown in Figure 3a, but the 
zonal axis has been expanded to be compatible with subsequent examples. 


The characteristic "wavelength" associated with this and subsequent examples 


is half that estimated earlier for flow on a plateau, that is, all stable 
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lownstream wave-like motions have a characteristic wavelength equal to that 
fa stationary Rossby Wave. 

Figure 5 shows the normalized streamline for the second case. The 
next example, as shown in Figure 6, indicates the same type of instability 
previously found for a single step with rp>2. For the assumed ro = 0.5, 
his type of instability exists for all 0.5 pa between’/2 and 37/2, 

Figure 7 shows the undisturbed downstream flow pattern that results 
when the plateau is an integral number of wavelengths long. This is the 
only topographic condition that allows the downstream flow to proceed 


undisturbed. 


Iv-3. Other Discontinuous Bottom Profiles 
We next consider the case where the flow is over a valley such that 
Yo = 1.5. The stability criterion (22) shows that such a flow is stable 
for all values of a. Figure 8 indicates the flow regime with pa = 5773. 
The flow for other choices of pa is similar but with different amplitudes. 
Positioning the step at integral multiples of 2773 gives an undisturbed 
downstream flow. The wavelengths of the disturbances over the valley are 
of course two-thirds the wavelength of the downstream waves. 

Returning now to the general solution, we recall] that the final 
depth was chosen equal to the initial depth solely for convenience. We 
may rewrite the solution (21) to include a fina] depth h3 = constant, 


where r3 * h3/H; 


Y 3(x,y) mo 4 1 + (1 - ro)sin ropa sin r 3p(x - a) 


(2h ) 
- f(a - ro) cos ropa + (ro/r, - 1) 2 cos r,p(x - a) 
2) cos ropa + (rg/ry = )} EB cos ryplx - a 
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For r3<r2<1, the flow is over two successive shoaling steps, and 


conversely, for r32r 21, the flow is over two successively deepening 


2 
steps. Intermediate alternatives represent modifications of the earlier 


examples. 


Iv-4,. The Bottom Profile is Smoothly Sloping 

Over any region of constant depth, the wavelength of the 
disturbance is fixed by the upstream velocity, the relative depth, and X. 
In order to investigate the possible influence of the mode of transition 
from one region of constant depth to another on the amplitude of the 
disturbance, it will be instructive to consider flow over a gradually 
shoaling slope. Returning to the general equations (10) through (16), 

we define the depth in the transition region as: h(x) = Fe) for XZ Xo. 
We are free to choose any numerical value of X,) that gives a desired slope. 
For simplicity we choose px, = 1, which gives an x, on the order of 100 km. 
This in turn describes a slope that halves the depth in about 100 km. 


Substitution then gives the transformed variables 





-e75 3 
D(s) = _ Ae ein oe. S,.%.0 ; i 
The solution of equations (12) and (15) is then 
Go = - 3UH(e®5 + cos s - sins). (25) 


In terms of the original variables, the flow over the slope has the forn 


Yo(x,y) 2 - SUHy +/x, + cos(1n x/x,) - sin(ln 1 %0)) (26) 
Bir px, yf 1, the solution is 


we 


| : PX epe(px,in x/x,) 4, 
Yo(x,y) Ts Poe fa es ad sin(px,1n x/x,) 5 (27) 








} 












; match solution (26) to flow on a flat plateau. We choose a 
of 0.4, which requires that we terminate solution (26) at 


' 2.5 X,- The matching conditions are then 


o(2.5x,) * £3(2.5x,) = Ag *°-1.15 UH 


o(2.5x,) = £3(2.5xQ) = r3pBz * -0.22 pUH; B; = -0.55 UH. 


ne complete solution is 


(x,y) = y $3(x) = - UHy [2.5 - 1.35 cos r3p(x-2.5x,) + O255 gam P3p(x-2.5x,)} 
(31) 


Figure 9 shows equation (31), again as a normalized streamline. 
Plotted on the same figure is equation (17), with r= 0.4. It is apparent 
that a transition slope of the above dimensions causes a phase shift but 
Joes not significantly change the amplitude of the deflection, which in this 
example is only 3 percent less for the flow including a transition slope 


than for flow without the slope. 
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CHAPTER V 


DISCUSSION 





V-1. Physical Interpretation of Results 
Sawyer (1959) classified the motions associated with topographic 


effects as follows: 

(a) Small scale turbulent motion mainly near the ground or ocean 
bottom. 

(bd) Gravity oscillations which may extend upward to great heights. 

(c) Oscillations about geostrophic equilibrium. 

(d) Quasi-geostrophic motions. 

The frictional effects denoted by (a) remain largely an unknown 

factor in large scale oceanic flows. Gtrtler (191) utilized the boundary 

layer principles developed by Prandtl to define an equivalent topography 

as that relief having the same effect on a frictionless current as does the 

actual topography on the real frictional current. He then treated the flow 

in the friction layer as a secondary current flowing normal to the main 

stream. A detailed analysis by GUrtler showed that both the magnitude of 

the secondary current and the deviation of the equivalent, or effective, 

topography from reality are usually so slight as to be negligible. Stommel 

(1957) argued convincingly for the treatment of Ekman convergence and 

divergence at both the surface and bottom of the sea in the same manner as 

precipitation and runoff are treated, namely, as external parameters 

influencing the internal flow regime. In any case, the topographic features 


under discussion are taken to be large enough that they do not directly 
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influence the frictional processes. There must be a point, however, at 
which the topographic disturbance of an originally frictionless flow may 
become large enough to require the consideration of friction in the 
subsequent flow. This is likely to be the case for the unstable flow 
patterns shown in Figures 3b and 6. 
By the assumption of homogeneity gravity oscillations of other than 
the free surface have been precluded, although it is recognized that 
internal oscillations may be an important, or even a dominant, result of 
topographic perturbations. If we adopt Sawyer's (1959) definition of 
gravity oscillations as those with pericds much less than a half pendulum 
day, we readily see that they may be neglected for our steady barotropic 
model. The shallow water waves associated with the depths of the ocean 
have phase velocities on the order of 100 m/sec, certainly faster than any 
current. Conversely, no deep water gravity waves exist with a phase 
velocity of 10 cm/sec (a reasonable choice for a current velocity). Thus 
it seems reasonable not to expect any stationary, surface, gravity waves 
to be associated with topographic effects on deep ocean currents. 

The classifications (c) and (d) appear to merge in a barotropic 
model. It is essential to note that the term "quasi-geostrophic" is a 
result of the utilization of perturbation techniques, and may be misleading. 
To illustrate this, Gambo (1957) demonstrated that the quasi-geostrophic 
_ approximation is poor for surface slopes of greater than 1/1000, or about 
5 minutes in angular measure. This is certainly less slope, by an order of 
magnitude or more, than is frequently found for both subaerial and submarine 
topographies. Thus it is believed that an examination of the developments 


given in the previous three chapters may give some insight into those 












topographic effects that lie without the region of applicability of 
linearized perturbation techniques. 

Returning to the fundamental equation for the conservation of 
potential vorticity, equation (lL), the qualitative effects of shoaling 
water (or deepening water) on an initially uniform zonal flow may be 
readily seen. For flow to the east, shoaling conditions immediately 
require a decrease in relative vorticity, or the establishment of 
anticyclonic flow. This turns the streamline towards lower latitudes 
where the Coriolis parameter is also decreasing. Depending on the change 
of depth, a point will be reached where the decrease in the Coriolis 
parameter will have balanced the change in depth, and subsequent motion 
towards the equator will cause the flow to become cyclonic thus returning 
the streamline to its original latitude. The converse is true for eastward 
flowing currents passing over deepening relief, except that there is no 
theoretical limit on the magnitude of the depth change, and hence equilibrium 
may not be reached. 

If the flow is from the east however, equation (l) shows that a 
wave-like motion is not expected, because there is no "restoring force." 
Shoaling conditions again require the development of anticyclonic flow, 
which in this case deflects the streamline to higher latitudes, and therefore 
higher values of the Coriolis parameter. This in turn requires a still 
further decrease in relative vorticity. The general solution (16) is 
applicable to flow in either direction. For eastward flows, s is real, 
giving a streamfunction periodic in s and therefore in x. For westward 
flows, U is negative, giving an imaginary value for p. This then transforms 
the trigonometric terms in solution (16) to give an exponential dependence 


of the streamfunction on x. 
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This qualitative analysis of tne vorticity equation is equally 
applicable to the Northern and Southern Hemispheres if the y axis is taken 
poleward, the x axis eastward, and the sign of the relative vorticity is 
defined as agreeing with the sign of the local vertical component of the 


earth's rotation. 


V-2. Comparison with Previous Investigations 
It is interesting to note that the results of this paper may be 


qualitatively compared with those of Fofonoff (1962) for the steady inertial 
circulation of a rectangular ocean of uniform depth and barotropic 
properties. Fofonoff obtained solutions to equation (0) of this paper 
(without the depth variation) subject to the condition that the boundary 
be a streamline. His solution for positive U is also periodic in x, but 
he chose to interpret it as a time-dependent solution superimposed on a 
Steady zonal current. 

The occurrence of critical flow patterns was discussed by Bolin (1950). 
He utilized the perturbation equations to investigate the flow of a very 
narrow zonal current over a Single mountain range. A critical velocity was 
found for which the flow would return to the source region. That velocity 
was given as the ratio of the local value of the Coriolis parameter to the 
total neight of the atmosphere, multiplied by the cross sectional area of 


the mountain range. Although the relevant choice of parameters was at or 





beyond the limit of the perturbation equations, Bolin demonstrated that 

such a critical regime might exist for features on the order of the Rocky 
Mountains. It is apparent that the streamlines depicted for critical flow 

in Chapter IV do not return to the source region because of the specification 


of an infinitely wide upstream zona] flow. 
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V-3. Numerical Evaluation and Comparison 
(a) iessabor ic Example 

Queney (1949), in his treatment of a stably stratified atmosphere 
of infinite height, gave the perturbations associated with flow at 10 m/sec 
over a plateau of height 0.5 km and width 300 km. Although the lack of a 
finite height of the air column in Queney's analysis hampers comparison, 
his results may be qualitatively related to those of Chapter IV by accepting 
Bolin's (1950) statement on the equivalent barotropic atmosphere for the 
case of stable stratification. Queney found only vertical perturbations 
of the type (b) above. He stated that lateral deflections are negligible. 
Using 8 km as the undisturbed height of the atmosphere, and the above 
dimensions of the plateau, the flow over the plateau is described by equation 
(17) as 


-¥o(x,y) = 8.5 x 10%y (1 - 0.06 cos 1.3 x 1075x). 


Tne coefficient of the cosine term in the above equation would seem 
te incleats a small deflection, but recalling that the y axis originates 
at the equator, this coefficient would give a deflection of about 6 percent 
of the original latitude were the plateau broad enough. However, the 
argunent of the cosine term varies through only 22 degrees for a breadth 
of 500 km, giving a maximum deflection at the leeward edge of the plateau 
of about 0.4 percent of the original latitude, or about 10 km at 30°N. 

From equation (21), downwind of the plateau 

-Y; = 8.0 x 108 E + 0.00) cos p(x - 500 km) + 0.02 sin p(x - 500 Kn) 
From this equation it is apparent that the amplitude of the downstream 
deflection is approximately 2 percent of the original latitude, or about 
79 km at 30°N. This amplitude must be compared with a representative 


wavelength of about 5000 km, and. therefore can be considered small. 
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Considering the above example, if the results of Chapter IV are to be 
applicable to atmosphcric flows, the scale of the relief must be larger 
than that considered by Queney. In fact, it was shown in Chapter IV that 
a more representative wavelength of a downstream oscillation is on the order 
of 10,000 km, and it is therefore clear that the only features that will 
give an oscillation of appreciable amplitude are those having a character- 
istic zonal dimension that is a significant fraction of the wavelength of 


the free wave. 


(b) Gulf Stream Meanders 

If we now attempt to compare the results of Chapters II - IV with 
observed oceanic currents, we are met by a marked lack of data for the 
barotropic mode. Rossby (1938) pointed out that, although there will 
always be a tendency towards the redistribution of mass in the oceans that 
will reduce horizontal pressure gradients with depth, such adjustments are 
not immediate and changes in the distribution of wind stress and atmospheric 
pressure must induce pressure gradients at the floor of the ocean. Such 
transient barotropic responses have not been included in this paper's 
steady state model. Recalling that the existence of barotropic flow was 
postulated earlier, the presence of oxygenated water in the deep ocean 
basins is accepted as primary evidence that a net flow extends to the depths 
of the oceans. 

Swallow (1957) found 0.8 cm/sec as the average velocity at about 
3000 m depth in the Northeast Atlantic. Volkmann (1962) measured deep 
currents on the order of 10 cm/sec at 3000 m just inshore of the Gulf Stream 
south of Cape Cod. Defant (1961) summarized earlier calculations of 
relative currents in the deep layers of the Atlantic. These calculations 


show velocities varying from 0 to 20 cm/sec over a range of depths from 
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2000 m to near the bottom. The wavelengths of Stationary Rossby waves 
associated with these velocities are on the order of a few hundred kilometers, 
certainly well within the existing scales of bathymetric relief in the deep 
oceans. It would therefore seem clear that the effects of topography on 
deep currents in the ocean should be appreciable over a much broader range 
of scales than those found to be inportant in atmospheric flows. | 

The solution (16) may be evaluated, if only numerically, for any 
relief of sufficient interest. In view of the rather restrictive assumption 
of an initially uniform zonal flow, and the lack of data on in situ currents, 
such an evaluation should await further refinement of the theory. One 
problem that could provide sufficient interest for future evaluation is 
that of the Gulf Stream; specifically, the problem of the meanders that 
appear after the Stream departs the Blake Plateau. 

Stommel (1960) discussed the meanders of the Gulf Stream at length. 
He utilized a perturbation technique to argue tentatively that the instability 
of such meanders may be associated with a critical internal Froude number. 
As another possibility he showed that the potential vorticity of the Stream 
is constant over a cross section, and thereby proceeded to develop a 
linearized theory of stable meanders. Stommel also mentioned the 
possibility that sea mounts may influence the development of the meanders, 
but did not develop a mechanism for such an influence. 

Such meanders are characteristically about 200 km in wavelengtk. 
Although any connection between the developments in Chapters III and IV 
and the real highly baroclinic Gulf Stream are at this stage tenuous at 
best, it is interesting to calculate the barotropic flow that would be 
associated with meanders of the dimensions observed. For L * 200 km, 


and again taking # = 2°10713sec™ tem}, the numerical relationship 
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ro * O.5 Uz( cm/sec) is obtained. Observing a rather rapid increase in 
depth from 1000 m to 3000 m in the region where the Stream leaves the Blake 
Plateau, ro is taken as roughly 3, in keeping with the introductory remarks 
concerning the utilization of realistic bathymetry. The instability 
associated with this choice of ry will be neglected initially, but will 

be discussed in a subsequent section. Such a change of depth would produce 
meanders of the observed wavelength in a barotropic current flowing at 

18 cm/sec. Volkmann (1962), although primarily interested in the apparent 
counter-current inshore of the Gulf Stream, showed (Volkmann, Fig. 3) the 
18 cm/sec isopleth at a depth of 2200 m with flow in the same direction 

as the Stream. The location is well removed from Cape Hatteras, being at 
approximately 36°N - 68°W, but may indicate strong barotropic currents in 
the Gulf Stream region. 

A more moderate change of depth would require a smaller velocity. 

That such physically ee are required despite the rather 
serious violation of the initial assumption of a uniform zonal flow indicates 
that such an assumption may not be too essential. 

As has been noted above, there is at present no convincing physical 
argument that allows the supposition of an equivalent barotropic Gulf Strean. 
Only if there is a degree of coupling between the surface Stream and the flow 
at depth can the above results be interpreted as indicating a possible cause 
and effect relationship. It seems reasonable to infer, however, that if a 
barotropic flow exists independent of the Gulf Stream, and it is subjected 
to the above described change in depth, the surface current may be influenced 


by the flow of the deeper layer. 
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(c) South Atlantic 

we now turn our attention to the South Atlantic and attempt to 
compare the results of Chapter IV with some of the observed features of the 
dynamic topography considered by Shtokman (1947) (see Figure 2). Retaining 
the Lagrangian viewpoint, Figure 2 shows an eastward flow departing the 
broad continental shelf in the vicinity of 40°S. The subsequent southerly 
excursion of the transport isopleths is approximately centered over the 
rather broad Argentine Basin (depth>5UL00 m). This is in qualitative 
agreement with the results of Chapter IV, which give poleward deflection 
over deepening relief. 

If, as for the Gulf Stream, a relative depth of 3 is assumed (again 
neglecting for the moment the instability implied by equation (18)), and 
observing that the wavelength along 40°S is about 8 degrees of latitude 
(about four times that of the Gulf Stream meanders), the velocity required 
is found to be 16 times that required for the barotropic component of ine 
Gulf Stream, or about 3 m/sec. This certainly bears no relation to 
physical reality. 

Similarly, downstream of the Argentine Basin the wavelengths of the 
dynamic contours remain on the order of 10 degrees of latitude, which for 
deepening relief require unrealistic velocities. Shoaling conditions however 
(say ra3), could give realistic velocities for such wavelengtns. The 
complexity of the bathymetry and of the upstream flow pattern precludes a 
realistic evaluation utilizing the model of Chapter III - IV. It should be 
noted that the apparent discrepancy cited by Shtokman (1947) between models 

in wnich the direction of the deflection depends only on the hemisphere, and 
, the dynamic topography shown in Figure 2, is subject to re-interpretation. 


Considering the region over the Mid-Atlantic Ridge at )0°S - 10°W, Shtokman 
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felt that the general southerly set of the current over gradually shoaling 
water contradicts Ekman and Sverdrup. Defant (1961) however noted that the 
sharp anti-cyclonic curvature of the isopleths of dynamic height approxi- 
mately over the top of the Ridge substantiates the Coriolis parameter 
models. The solutions obtained in the previous chapters of this paper of 
course require equatorward deflection over shoaling water for upstream 
zonal flow, but do not explicitly treat the probable interactions between 


the effects of successive bathymetric features. 


(d) Antarctic Circumpolar Current 

The Antarctic Circumpolar Current (ACC) affords the classic example 
of topographic effects on zonal currents (see Figure 1). In general the 
relative (baroclinic) currents are observed to be deflected towards the 
equator over shoaling water and towards the pole over deepening water, 
an observation which qualitatively agrees with our results. The sharp 
northward perturbation of the isopleths over the New Zealand Plateau will 
be considered first. The apparent lack of any downstream effects indicates 
that the Plateau should be about one wave length long, if the results of 
Chapter IV are to be applicable. The relative depth associated with the 
Plateau is about 0.1 and a characteristic width at 5095 is 780 km. These 
figures require a barotropic velocity of about S cm/sec, a reasonable value. 

The flow of the ACC over the Scotia Ridge and South Sandwich Trench 
system affords an even more striking example of topographic effects, although 
the constrictive effects of Drake Passage tend to complicate the upstream 
flow pattern. The complexity of the bathymetry in this region, including 
the marked cross-stream variations in depth associated with the South 


Sandwich and Orkney Islands and South Georgia precludes all but the roughest 
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estimates of relative depth and feature breadth. Such an estimate however 
does yield a physically reasonable velocity on the order of 1C cm/sec. 
Between 60°© and 90°E Figure 8 shows that the ACC is deflected 
towards the equator as it passes over the Kerguelen - Gaussberg Ridge. 
Again utilizing the final investigation of Chapter IV to validate the 
approximation of such a smooth feature by a plateau, the distance between 
the 3600 m isobaths along 55°S is chosen as the representative width; 

about 550 km. Representative depths of 4.5 km upstream and 3.6 km over the 
feature with a relative depth of 0.8 are chosen. These values of the 
parameters require a barotropic velocity of 9.8 cm/sec for undisturbed 


downstream flow. 


V-4. Comparison and Discussion of Amplitudes 
It would be advantageous to be able to use the amplitudes of the 
observed disturbances aS independent cnecks of the agreement between theory 
end observation. Agreement by this means can not be obtained for the 
examples cited above. As has been noted previously, the relative depth of 
3 assumed for the Gulf Stream and for the Argentine Basin examples pive in 
éach case unstable flow conditions for tne frictionless model. The 
amplitude of the disturbance over the New Zealand Plateau should be 
approximately hO percent of the original latitude, or about 20 degrees, 
rather tran the observed 5 - 10 degrees. The calculated amplitude over 
the Kerguelen - Gaussberg Ridge is about 17 degrees as opposed to the 
observed amplitude of 7 - @ degrees. 

There are several possible explanations for the discrepancies between 
the amplitudes of the observed and calculated disturbances. It was shown 
at the end of Chapter IV that the mode of transition in depth from one region 


to another has little effect on the amplitude of the disturbances. Therefore 





34 


the step-shaped features considered are not the primary source of the 
disagreement. The departure of the real current from tne assumption of 
meridional uniformity is certainly appreciable in the case of the Gulf 
Stream, but would seem to be minor in the case of the Antarctic Circumpolar 
Current. 

Of the two remaining major factors, friction and baroclinicity, it 
seems probable that the latter is of fundamental importance in determining 
the amplitude of the —— disturbance. Saint-Guily (1962) noted 
that this difference in amplitude is likely due ".. . to the existence of 
a vertical density gradient. Such a gradient apparently has the effect of 
intensifying the currents near the surface, therefore of reducing the 
equivalent depth H, ... . , and the amplitude of the deflection." Our 
introductory postulation of the existence of a barotropic mode of oceanic 
flow included no a priori expectation that such a flow over physically 
realistic topographic features would directly correspond to the observed 
relative dynamic topography. It might be expected that an equivalent mode 
of circulation with uniform velocities extending to great depths would also 
require the introduction of the concept of an equivalent topographic feature. 
Unfortunately there is as yet no independent procedure available by which the 
relief of such a feature may be inferred, 

In the brief examples given on previous pages of the chapter the 
observed wavelengths and the actual bathymetry have been utilized to calculate. 
barotropic velocities. This approach may be modified and both observed 
wavelengths and observed amplitudes utilized to calculate equivalent 
velocities and topographies. For the flow over the New Zealand Plateau 
shown in Figure 1, the amplitude is about 10 degrees of latitude and the 


wavelength is about 780 km. Substitution of these values into equation (18) 





op) 


yields an equivalent relative depth rj, of about 0.9 (as opposed to the 
physical change in depth of about 0.4). The resulting barotropic velocity 
is about 25 an/sec. 

A similar evaluation of the flow over the Kerguelen - Gaussberg Ridge 
yields an equivalent relative depth of 0.94 and a barotropic velocity of 
1h cm/sec. 

A characteristic amplitude of the Gulf Stream meanders is less 
easily defined, but as a rough approximation it may be taken to be of the 
same order as the wavelength, about 200 km (Stommel, 1960). This value 
gives a barotropic velocity of 2 cm/sec, and thé very smal] relative depth 
Ono. 03. 

The above consideration of an equivalent topography indicates the 
reason for the qualitative success of the perturbation theories. Such 
theories would be seemingly justified on the basis of the small effective 
changes of depth that were above found to be required to produce disturbances 
of the barotropic mode that correspond to the observed disturbances of the 
baroclinic component of the ocean currents. Such a distinction between 


actual and effective topography has not generally been made. 
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CHAPTER VI 
SUMMARY 


1. An exact solution has been obtained for the steady disturbance of 
an initially uniform, frictionless zonal] barotropic current by a meridionally 
uniform topographic feature. 

2. The resulting streamlines are zonally periodic for eastward flow 
and exponential for westward flow. 

3. In both the Northern and Southern Hemispheres eastward flow is 
deflected towards the equator over shoaling water and towards the poles over 
deepening water. 

lh. Topographic effects are not limited to the region of variable 
depth, but may exist downstream in the form of a standing Rossby Wave. 

S. For a finite decrease in depth between two regions, the manner 
of transition in bathymetry is unimportant and may be treated as 
discontinuous. 

6. The wavelengths of stationary Rossby waves associated with 
characteristic velocities of the deep currents in the ocean are such that 
the existing scales of topographic features should be highly effective in 
disturbing deep barotropic currents. 

7. Very small changes in depth in the barotropic model are sufficient 
to produce disturbances of the dimensions observed in the real ocean. 

8. Future study should consider the effects of the baroclinicity of 
the ocean, which should be properly related to both the actual topography 


~ 


and the vertical structure of the oceanic currents. 
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Figure 3b. Zonal flow over a deopening step. 
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